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Abstract 

We propose a program for finding the cyclicity of period annuli of quadratic 
systems with centers of genus one. As a first step, we classify all such systems 
and determine the essential one-parameter quadratic perturbations which pro- 
duce the maximal number of limit cycles. We compute the associated Poincare- 
Pontryagin-Melnikov functions whose zeros control the number of limit cycles. 
To illustrate our approach, we determine the cyclicity of the annuli of two 
particular reversible systems. 

0. Introduction 

As well known, there are four types of planar quadratic systems with a center: 1) 
Hamiltonian, 2) reversible, 3) generalized Lotka-Volterra, and 4) of codimension four. 
The first integral in the Hamiltonian case is a cubic polynomial in (x, y) whose generic 
level sets are elliptic curves. It is easy to observe that the generic level sets of the 
first integral of the codimension four case are elliptic curves too (see the end of Sec- 
tion 4). Our main efforts will be devoted to the remaining two cases. It has been 
recently proved by one of the authors (S.G.) that there are 18 classes of reversible 
centers, 6 classes of reversible Lotka-Volterra centers and 5 classes of generic (non- 
reversible) Lotka-Volterra centers whose phase portraits contain only elliptic curves 
(possibly, a finite number of them reducible) [5]. They are given by codimension-one 
or codimension-two algebraic sets in the space of all centers from the corresponding 
type, see Theorems 1 and 2 below. Throughout the paper, by "genus" we mean the 
genus of the compactified and normalized generic phase curves. An algebraic phase 
curve is generic if it does not contain a singular point of the vector field in its closure. 
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The centers whose (generic complexified) periodic orbits are elhptic curves will be 
called centers of genus one. We note that even a quadratic system can have a center 
of arbitrarily big genus. 

Once we know that the first integral of a given planar system with a center de- 
fines elliptic curves, we could raise the next question: how many limit cycles could be 
produced in the phase portrait under small quadratic (or even polynomial) pertur- 
bations? The purpose of the paper is to present a program for solving this problem. 
In what follows we restrict our attention to the limit cycles which bifurcate from 
open period annuli (we do not consider graphics). It turns out that, instead of 
multi-parameter perturbations, it is enough to consider suitable one-parameter small 
quadratic perturbations, see [9]. The one-parameter perturbations which produce 
the maximal number of limit cycles in the quadratic case, called essential perturba- 
tions, together with the corresponding generating functions of limit cycles (or also 
Poincare-Pontryagin-Melnikov functions) were determined in |14j . 

The first part of our program is to adapt the results of to our case. The result 
is a complete list of such essential perturbations of quadratic systems with centers of 
genus one, together with the corresponding generating functions. The list is presented 
in Section 3 (in the reversible case) and in Section 4 (in the Lotka-Volterra case). 

The second part of the program is to study the zeros of the generating functions 
I{t) found in Sections 3 and 4. The fact that each J(t) = is a complete 

elliptic integral from a rational one-form uj., over the level sets {H = t} C which are 
elliptic curves, allows one to apply all related facts from algebraic geometry in order to 
estimate the number of zeros of I{t) and thus to set up some upper (or lower) bounds 
on the number of bifurcating limit cycles in the system, see e.g. [201 [IHl [HI E E]- 
This part of our program is illustrated in Section 5, where we study the cyclicity of 
the period annuli for two types of quadratic reversible systems with centers of genus 
one. 

The paper is organized as follows. In Section 1 we determine all reversible centers 
with phase portrait formed by elliptic curves. In Section 2 we discuss the same ques- 
tion for the Lotka-Volterra centers (reversible or not). These results were previously 
proved in [5]. For convenience of the reader we present here almost self-contained 
proofs adapted for the purposes of the present paper. In Sections 3 and 4 we de- 
termine for each of the cases the corresponding generating function, the complete 
elliptic integral I{t) which is the leading term in the expansion of the first return 
mapping, respectively for the reversible and Lotka-Volterra cases. We also present 
there (as Conjectures 1 and 2) the expected exact upper bounds for the number of 
zeros of all generating functions. In Section 5 we use the geometric properties of the 
elliptic fibration determined by the first integral of the quadratic system, in order to 
determine, for two of the cases, the number of the zeros of I{t). From this we deduce 
an exact result: the cyclicity of the annuli under consideration is two (Theorem 3). 
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1. The reversible case 

The general first integral of a reversible system with a center at the origin 

z = —iz + az'^ + 2\z\'^ + bz'^ , a, 6 G M, z = x + iy (1.1) 

is given by 

H{x,y)=X\ly^ + AX^ + BX + C) (1.2) 

where X = l+2{a—b)x and X,A,B,C are explicit rational functions of the parameters 
a, b (see formula (1.7) below). Moreover, one has A 7^ 0, —1, —2. Along with the elliptic 
Hamiltonian H = + — (corresponding to a = 6 = —1 in (1.1)), formula 
(1.2) contains all the cases for which {H = t} is (possibly) an algebraic curve of genus 
1. Below, we are going to determine all these cases. 

We begin with the observation that it is enough to consider only the case when 
A < —1. Indeed, if A > —1, by the bi-rational change {X,y) = (l/Xi,Fi/Xi) the 
function (1) reduces to Hi = X^^-^{^Y^ + CXf + BXi + A) with -2 - A < -1. For 
this reason below we will consider the case when 

P 

A = — with p, g G N, p > q, p 7^ 2g, gcd(p, q) = I 

q 

where gcd{p, q) denotes the greatest common divisor of p and q. 

1. Assume first that A 7^ 0, C 7^ 0. Then {H = t} after taking X = (this is an 
isomorphic map [5], Lemma 1) reads 

ly-i = -Ax^'^ - Bx'^ - C + txP. (1.3) 

For arbitrary t, the curve (1.3) has a genus 1 if and only if the polynomial on the 
right hand side has a degree 3 or 4. That is, either p = 3 and g = l,2orp = 4 and 
q=l. 

Consider now the degenerate cases when either A = or C = (let us note that 
neither two of the coefficients A, B, C can vanish simultaneously). 

2. Assume that A = 0. Then (1.3) is of genus 1 if and only if either p = 3, g = 1, 2 or 
p = 4, g = 1, 3. Thus the unique new case compared to the case A 7^ is (p, q) = (4, 3). 

3. Let C = 0. If 2g > p > g, by a bi-rational transformation {x,y) = {l/xi,yi/x\) 
one can transform (1.3) into 

1^2 = _ Bx'i + tx^'-P (1.4) 

(here and below, we shall omit the subscript 1). As q > 2q — p, the genus is 1 if and 
only if either g = 3, p = 4, 5 or g = 4, p = 5, 7. If p > 2g and g is even, then the 
rational change y = yix'^/'^ in (1.3) yields 

^.y"^ = -Ax'i - B + txP-'i. (1.5) 
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As p — q > q, (1.5) is of genus 1 if and only if {p, q) = (5, 2). Finally, if p > 2q and q 
is odd, the rational change y = in (1.3) yields 



ly2 = -Ax'i^^ -Bx + txP-'^+\ 



;i.6) 



As above, (1.6) is of genus one if and only if p — g + 1 equals 3 or 4, which is possible 
when either {p,q) = (3,1) or {p,q) = (4,1), cases that already have been obtained 
when we assumed that A ^ 0. 

Thus we have obtained the complete list of cases for which {H = t}, H given by (1.2) 
and t arbitrary, is a curve of genus one (the right column contains the cases with 
A > -1): 



i)H = X-3(|j/2 + AX^ + BX + C) 
iii) H = X-3/2(i^2 ^ j^x^ ^BX + C) 
v) H = X-\\y^ + AX'^ + 5X + C) 
vii) E = X-4/3(i|/2 ^BX + C) 
ix) H = X-4/3(i^2 ^ ^^2 ^ BX) 

Xi) H = X"5/3(1^2 ^ ^ 

xiii) H = X-5/4(ly2 + AX^ + BX) 
xv) = X-'^'Wy^ + + BX) 

XVii) = X-5/2(1^2 ^ ^j^2 ^ ^j^) 



We should mention that in cases (iii) and 
C = 0, the curve has a genus zero, see below 

Let us now recall the exact formula of (1.2 



the standard elliptic case) 
ii) H = X(|t/2 + CX^ + BX + A) 
w) H = X-i/2(iy2 + cx^ + BX + A) 
vi) H = X2(ly2 + CX^ + BX + A) 
viii) H = X~2/3(i^2 + QX^ + BX) 
x) = X-2/3(ly2 + 5X + A) 
xii) = X-V3(i^2 ^BX + A) 
xiv) = X-3/4(i^2 ^BX + A) 
xvi) i7 = X-i/4(|y2 + 5X + A) 
xviii) if = Xi/2(|y2 + 5X + A). 



iv) above it was assumed that C ^ 0. If 



from 



a + b+2 

H{X,y) = X-^^ 



— + 

2 



1 



6)2 Va-36-2 



X^ + 2- 



-X + ^ 



36 + 2 



6 + 2 



(1.7) 

This formula holds outside the lines a = 6, a + 6 + 2 = 0, 6 = — l,a — 36 — 2 = (note 
the last three cases correspond to A = 0, —1, —2 in (1.2)). Except for the three points 
(a, 6) = (—1, —1), (±2, 0), on these lines the first integrals contain exponents [H] and 
hence their level sets are not algebraic curves. For (a, 6) = (±2, 0) the level sets are 
conic ovals. The lines a + 6 + 2 = — A(a — 6) in the (a, 6)-plane, A G M, together with 
a = b, form the bundle of straight lines through the point (—1,-1) which corresponds 
to the standard elliptic case. Therefore using the above results, we obtain: 



Theorem 1. The phase curves of {1.1) are algebraic curves of genus one if and only 
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if one of the conditions holds: 



(rl) a = 26 + 1 (r2) a = — 1 (the reversible Hamiltonian case) 

(r3) a = 56 + 4 (6^-3) (r4) a = -36 - 4 (6 ^ -3) 

(r5) a = §6 + I (r6) a = |6 - | 

(r7)(a,6) = (|,-i) (r8) (a, 6) = (-1, -i) 

(r9) (a, 6) = (-8, -2) (rlO) (a, 6) = (4, -2) 

(rll) (a, 6) = (-17, -5) (rl2) (a, 6) = (7, -5) 

(rl3) (a, 6) = (-7, -|) (rl4) (a, 6) = (^, -|) 

(rl5) (a, 6) = (-23, -7) (rl6) (a, 6) = (9, -7) 

(rl7) (a, 6) = (13, 5) (rlS) (a, 6) = (-3, 5). 

For completeness, below we add the list of all cases in (1.1) for which the ovals {H = t} 
are conic curves (ellipses). These are 

{rl9) H = X-^/^{^y^ + AX^ + BX) when (a, 6) = (-11, -3) 
(r20) H = X-i/2(i^2 ^BX + A) when (a, 6) = (5, -3) 
{t21) H = X-\^y^ + BX + C) when (a, 6) = (2, 0) 
(r22) H =ly^ + CX^ + BX when (a, 6) = (-2, 0). 

Cases (r21) and (r22) (not included in (1.2) and (1.7)) are taken from the full hst of 
first integrals of (1.1), see e.g. [H]. The cases (rl9) and (r20) are obtained from (1.2) 
in the same way as above. Note that (r20) and (r21) are the isochronous centers 
and 1S2, respectively. By the way, the isochronous center 1S4 corresponds to 6 = 1 in 
(r5). Perturbations of the quadratic isochronous centers have been studied in [3]. 

2. The Lotka-Volterra case 

A (generalized) Lotka-Volterra system with a center at the origin has in complex 
coordinates the form 

z = -iz + Az"^ + Bz^ , z,A,BeC. (2.1) 

Apart of the classical Lotka-Volterra system, the generalized one splits into real and 
complex cases. In appropriate coordinates, the general first integral of the Lotka- 
Volterra system in the classical real case is 

H{x,y) = x^yf'il -X -y), A,/iGM, A;u(A + ^ + 1) ^ 0. (2.2) 

The general first integral in the complex case is 

H{x,y) = {x^ + yy{l-2x)e-^^^"'"'''^y/''\ A,/ieM, A < (2.3) 

When /i = in the complex case and when (A — fi){X — — 1) = in the real 
case, the corresponding system becomes reversible. That is, after a suitable affine 
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change of the variables, the initial system with a first integral (2.2) or (2.3) will take 
in complex coordinates one of the normal forms 



z ^ -iz + z"^ + bz"^ , beR] z^-iz + z'^. (2.4) 

This is possible provided that the coefficients in (2.1) satisfy A^B G M. Our main 
result in this section is the following. 

Theorem 2. The phase curves of {2.1) are algebraic curves of genus one if and only 
if one of the conditions holds: 

(rlvl) ^4 = (the Hamiltonian triangle) 

(rlv2) 2AB - ^ (Ivl) AB + {1 ± 2i)A^ = 

(rlv3) AB - = (lv2) - (101 ± 2%i)A^ = 

(rlv4) bAB - 3^2 = (lv3) 289AS - (151 ± m)A^ = 

(rlv5) bAB - ^2 = (lv4) 1681^45 - (783 ± 60^2^)^^ = 

(rlv6) ?,AB + ^2 = (lv5) 841A5 - (349 ± 12^)^^ = 

The above statement is a consequence of the following two propositions which will be 
proved together: 

Proposition 1. The phase curves of the reversible Lotka-Volterra system (2.1), 
A^B e M, are algebraic curves of genus one if and only if its first integral is affine 
equivalent to one of the normal forms 

(rlvl) H = xy{l — x — y) (the Hamiltonian triangle) 

(rlv2) H ^ x~^y{l - x - y) 

(rlv3) H = x^y{l — x — y) 

(rlv4) H = x~^|/(l — X — y) 

(rlv5) iJ = (x^ + _ 2x) 

{r\v6) H = {x^ + y^)-\l-2x). 



Proposition 2. The phase curves of the generic Lotka-Volterra system (2.1), A'^B ^ 
M., are algebraic curves of genus one if and only if the first integral is affine equivalent 
to one of the normal forms 



(Ivl) H = 


2 

X 


1/3(1- 


■x-y) 


(lv2) H = 


x~ 


-V(i 


-x-y) 


(lv3) H = 


x~ 


-6^3(1 


-x-y) 


(lv4) H = 


X' 


-V(i 


-x-y) 


(lv5) H = 


X' 


-3yi(l 


-x-y) 
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Proof of Propositions 1 and 2. Under the conditions in (2.3), tlie critical point 
(dJ+T'IaTt) ^ center. The origin is a focus for /x 7^ and a center elsewhere. 
Clearly, the phase curves defined by (2.3) could be elliptic only if /i = 0. Note that if 
fi = 0, the origin is a center for all A 7^ 0, but this center is reversible and the system 
can be transformed into the normal form (1.1). Here we study the Lotka-Volterra 
center outside the origin, existing for A < — |. 

Under the conditions in (2.2) (frankly, one should take H = xy{l—x—y)\x\^~^\y\^~^ 
there, but modules will be omitted thoroughly), there is a unique critical point 

(^ a+a1+i ' A+/I+1 ) outside the invariant straight lines x = 0, y = 0, x + y = 1, which is 
a center if and only if A/i(A + fi + l) > 0. In fact, there are two topologically different 
configurations having a center, the first one is obtained for A > 0, /i > and the 
second one corresponds to the parameters outside the first quadrant. In this latter 
case, one can take without loss of generality A < 0, ;U < 0. Indeed, if e.g. A < < /i, 
then applying an affine change x = 1 — X — Y , y = Y, we reduce the first integral 
H{x,y) = t in (2.2) to Hi{X,Y) = ti where 

Hi{X,Y) = X^/^Y^'/^{l- X -Y), ti = t^'^ (2.5) 

One can proceed similarly with the other case yU < < A. In the same way, when A 
and /i are positive, we can reduce their values to A,/i G (0, 1]. Indeed, if e.g. A > 1 
and A > /i, the same change as above transforms (2.2) into (2.5), with both degrees 
in (2.5) within (0, 1]. We proceed similarly if yU > 1 and /i > A. 

When A/i(A + yU + 1) = 0, the first integral in (2.2) should be replaced by another 
one containing logarithmic or exponential terms [HI [211 123] and therefore its level 
sets are not algebraic curves. In this way we have reduced our consideration to the 
following two cases 

(I) < A < 1, < yU < 1, 

(II) A<0, yU<0, A + yU+l>0. 

In the reversible Lotka-Volterra cases, one can find all curves with a genus one by 
simply using the results from the previous section. Take the real reversible case and 
assume for definiteness that A = /i. The remaining possibilities are reduced to this one 
by the same change which we used to obtain (2.5) from (2.2). Let A = yU G (—1/2, 1]. 
A further substitution in (2.2) x = |(1 — X) +Y,y = |(1 — X) — Y transforms the 
curve H{x, y) = t to Hi{X, Y) = ti where 

H,iX, Y) = XVA[_ir2 + _ 1)2]^ ^ (2.6) 

As 1/A G {—00,— 2) U [1,00) and (2.6) takes the form already studied in Section 1 
above, we conclude immediately that Hi, and hence H, is of genus one if and only 
if 1/A G {-4,-3,1,2}. Namely, for A = yU G {1, i, -|} in (2.2). In the same 
way, taking in the complex case (2.3) fi = and using the substitution x = |(1 — X), 
y = Y, we reduce (2.3) to (2.6) (with a sign + in front of |). As 1/A G (—2,0) 
now, one obtains a curve of genus one if and only if 1/A G { — |}- That is, for 
A G { — I, —2} and yU = in (2.3). Thus we have obtained the complete list of cases 
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with genus one in the reversible Lotka-Volterra system (the first three entries in the 
right column contain the cases with /i = 1; the cases with A = 1 are obtained by a 
rotation of the variables {x,y) {y,x)) 



H = 


xy{l — X - 


-y) 










H = 


x^y^{l — 


x-y) 




H = 


x'^y{l — X - 


y) 


H = 


(1 


— X — 


y) 


H = 


a;~'^?/(l — X 


-y) 


H = 


(1 


— s — 


y) 


H = 


x~^y{l — X 


-y) 


H = 


(X2 + 


-1(1- 


2x) 


H = 


(x^ + ?/2)-2| 


[l-2x 



For a completeness, let us add to the above list the unique reversible Lotka-Volterra 
case having conic orbits, namely 

(rlvO) H = {x^ + - 2x). 

This is the quadratic isochronous center known as Si and it corresponds to 6 = in 
(2.4). This finishes the proof of Proposition 1. 

Now we come to the generic (non- reversible) real cases (I), (11) above, where 
(A - /i)(A - l)(/i - 1) ^ 0. Consider first (1). Let 

A = -, /i = -, p,q,r,seN, p<q, r < s, gcd{p,q) = l, gcd{r,s) = l. 
q s 

We replace our first integral H{x,y) = t with Hi{x,y) = ti where: 

Hi{x,y) = xP^y''^{l- X -yY^, p^ = Xn, qi = nri,ti = f\ n = qs/gcd{q, s). 

(2.7) 

The irreducible algebraic curve 

Ft = {(x, y) G : x^y'^{l — x — yY = t}, p, g, r G N, p < q < r, gcd{p, q,r) = 1 

(2.8) 

2 I 2 I 2 

is smooth if and only if t ^ A = {0, ( ^^^_^_^ ) ? }. Let Ff be the associated compact 
Riemann surface. Its (geometric) genus is one and the same for all t ^ A. We have 
(cf. [51E]): 

Proposition 3. The compact Riemann surface Tt, t ^ A, is of genus one if and only 
^fp=^,q = '^,r = 3. 

Proof. For reader's convenience, we repeat the proof from [5]. It is based on the 
Poincare-Hopf formula, namely deg{ijj) = 2g — 2, which we apply to the meromorphic 
one-form related to the first integral (2.8) 

dx dy 

x[q — qx — {q + r)y] y[p — py — {p + r)x] 

We recall that, in our context, deg{uj) = Xlfc'^fc where the summation is taken over 
all zeros or poles of u, and Uk are their orders. The above one-form has neither zeroes 



8 



nor poles in the affine chart outside the critical locus. Hence it suffices to consider it 
at infinity. There are 3 points at infinity: 5*1 = [1 : : 0], 5*2 = [0 : 1 : 0], 5*3 = [1 : 
— 1 : 0]. The local coordinates near can be chosen as follows. Write x = 1/u with 
— > 0. After this change of coordinates, equation (2.8) becomes y'^{l—u+yuY ~ u^'^'^ 
(we take t — (—1)^ for simplicity). Since w — > 0, we have y ^ (if y — > oo, we would 
have yu ^ —1 which corresponds to S3). Consequently, up to bi-analyticity, there 
are local coordinates {u,y) verifying — u^^^ . Let m = gcd{q,p + r). There are m 
different parameterizations, namely: 

which correspond to the m different local coordinates near the m smooth points given 
after blowing-up 5*1. In these local coordinates, 

—ud{l/u) qd^ 
^ ~ q-q/u-{q + r)y~ m^[q _ q^-^ _ (g + r)e^^^] ' 

We conclude that near the m points, cu has a zero of order — — 1. 

We perform the same study near the two other points at infinity. Near 5*2, by 
symmetry, we obtain n = gcd{p, q + r) points where the one-form u has a zero of 
order ^ — 1- Similarly, near S3, there are / = gcd{r,p + q) points where lj has a zero 
of order y — 1. 

As a result, we obtain the formula deg{uj) — p + q + r — n — m — I and therefore 
the curve (2.8) is elliptic if and only if 

p + q + r = n + m + l. (2-9) 

Now we have to resolve this Diophantine equation. Obviously, we have m < q, n < p 
and / < r, hence (2.9) is true if and only if 

gcd{q, p + r)^q, gcd{p, q + r)^p, gcd{r, p + q)^r. 

Take natural numbers a, /3, 7 such that q + r = pa, p + r = q(3, p + q = rj. The 
latter system has a nonzero solution 

p _ 7 + 1 g _ 7 + 1 
r a + r P + 1 

if and only if 

a(3'y = 2 + a + P + -f. 

By (2.8) wc have a > /9 > 7 and the unique solution satisfying this condition is a = 5, 
/3 = 2, 7 = 1 which leads to p/r = 1/3, q/r = 2/3. □ 

To consider the second case (II) , we note that a bi-rational change of the variables 

_ -X _ -Y 

^~ 1-X-Y' ^~ 1-X-Y 
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transforms (II) to already considered: Hi = X^Y^{1 - X - Y) = ti with 

A > 0, M > 0, where 

A=-^^ , M= ~^ , =ri/(A+M+i). (2.10) 

X + fi + 1 X + fi + V ^ ' 

Hence one can use the values of A, M already obtained above to calculate A and [i 
for case (II) through the formula 

-A -M 



A+M+r ^ A+M+l 

Thus we have completed the list of all cases with genus one in the generic Lotka- 
Volterra system (another 15 cases are obtained by a rotation of the variables): 



H = 




1 — X 


-y) 




H = 


3 1. 

X2y2 [ 


1 - 


- X - 


-y) 


H 


= x'^y^\ 


[1- 


■ X — 


y) 


H = 


1 

X &y~ 




X — 


y) 


H = 


x^^y"^! 


(1 


— X 


-y) 


H 


1 

= X2y 


-3(1 


— X 


-y) 


H = 


x'ly- 




X — 


y) 


H = 


x^^y^' 


(1 


— X 


-y) 


H 


1 

= X3y 




— X 


~y) 


H = 


1 

X 


4(1- 


X — 


y) 


H = 


x^^y"^! 


(1 


— X 


~y) 


H 


1 

= X2y 




— X 


-y) 


H = 


1 

X 


4(1- 


X — 


y) 


H = 


X y2 


(1 


— X 


-y) 


H 


2 

= X3y 




— X 


-y) 



Proposition 2 is completely proved. □ 

Proof of Theorem 2. The proof is a direct consequence of Propositions 1 and 2. 
More precisely, it follows by a straightforward calculation using the formulas from 
case (iv) on page 157 there. □ 



3. The generating function in the reversible case 

We are going to study small perturbations of the reversible system (1.1) in the cases 
when the first integral H defines a curve of genus one. Consider first the quadratic 
perturbations of system (1.1) rewritten in real coordinates: 

X = Hy/M + ef{x,y,e), 
y = -H^/M + eg{x,y,e) 

where M = X^^^ and f,g are quadratic polynomials with coefficients depending an- 
alytically on the small parameter e. Given a perturbation {f,g), the limit cycles in 
(3.1) are determined by the zeroes of the leading term I{t) in the expansion with 
respect to e of the displacement map. For this reason the integral I{t) is called some- 
times the generating function. As well known |T3], one can always choose a particular 
quadratic perturbation so that the corresponding I{t) would have the possible maxi- 
mum of zeroes within the whole class of quadratic perturbations. In the generic case 
I{t) is given by the complete elliptic integral [H], Theorems 2 and 3 



I{t) = / x'^ ^ {^i + + fisx ^)ydx (3.2) 

J5{t) 
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where S{t) is the oval contained in the level set H = t and fii G M. In the exceptional 
cases (rlO) and (r5) with 6 = 2 when system (1.1) belongs to the intersection with the 
codimension-four stratum Q4 of the center manifold (see |23j for details), I(t) takes 
another form and is not an Abelian integral [13] • In the standard elliptic case when 
a = b = —1 in (1.1) (we will denote it by (rO)), the integral has a specific form, too. 
Using (3.2) and the list of first integrals of genus one (O)-(xviii) above, we obtain the 
concrete form of I{t) for all the cases (after appropriate re-scaling of x, y, t, H, I) as 
follows: 



(rO) 


I[t) 


r 

~ Jh=o 




TT 

£1 = 


1^,2 
2^ - 


(rl) 


I{t) 


~ Ih=o 




T (+\ 


f 

Jh=q 




u 

n = 


2y + 


(r3) 


m 


~ Ih=o 




I{t) 


f 

Jh=o 




TT 

n = 


2^ + 


(r5) 


m 


~ Ih=o 


{lb) 


T(+\ 


r 

Jh=q 


[lb) 


l[t) 


_ r 
Jh=o 




H = 


2^ ^ 


(r7) 


l[t) 


r 

Jh=o 


(r8j 


T(+\ 


r 

Jh=o 




TT 

H = 


1 2 1 
2^ + 


(r9) 


m 


~ Ih=o 


(rlO) 


lit) 


~ Ih=o 




H = 


v + 

2^ ^ 


(rll) 


m 


~ Ih=o 


(rl2) 


i{t) 


~ Ih=o 




H = 


\y' + 


(rl3) 


I{t) 


~ Ih=o 


(rl4) 


I{t) 


~ Ih=o 




H = 


\y' + 



/iy=o(/^i + Ai2t + fJ^3x)ydx, 

x~^{fii + fi2X + fi3X^^)ydx, 



X "^(/ii + fi2X ^ + fi3x)ydx 



3~b 



_L ti^ I l-3b „2 
6(fe+l) 6+1-^ 2(6+1)-^ 



tx^^ 



X '(yUi + iJ,2x'^ + nsx '^)ydx, 

2 I ,, ^2^1 



-2 
6+3 



a; (yUi + fi2X + )ydx, 



b-l 

24(6+1) ' 4(6+1) 



+ 



X^ + 



36+1 
8(6+1) 



X — tx 



^x ^{jji + j^2X + jJ^sx ^)ydx, {b^2) 
^x~^[fii + /i2X + fi3X~^ + /i4(l — x~^) lnx]ydx, {h = 2) 



X ^{fii + fi2X ^ + fisx)ydx, 



2-b 



L tli^ I 1-26 „2 

4(6+1) 6+1-^ 2(6+1)-^ 



tX^ 



X 

J_ 

12 



1„3 



^(/ii + /i2X^ + /isx '^)ydx, 
^"^{fii + fi2X~^ + fi3X^)ydx, 
-tx\ tG(-i,0) 

+ fi2X^^ + fi3X^)ydx, 

g + ^x^ - tx'^, t e (|, 00) 

,a;(/ii + /i2X"^ + fisx^)ydx, 
x~^{^i + /i2a;^ + fi3X~^)ydx, 

00) 



X 



*) In 



^2' 



/ff=o(/^i + y"2a;"^ + fi3X^)ydx, 

^ x"^(/ii + /i2a;'^ + iJ,3X~^)ydx, 
y^ + ^ + lx^-tx\ te (|,oo) 
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(rl5) I{t) = /f^^oX^(/ii + /i2a;~'* + fi3X^)ydx, 
(rl6) I{t) = J^^QX~'^{j^i + j22X^ + j^3X~^)ydx, 
H =ly^ + j + j2^'^ -tx, te (|, oo) 

(rl7) I{t) = /^^Q x~^{fii + /i2x2 + ii3x~'^)ydx, 
(rl8) I{t) = Jjj^QX^^{fii + fi2X^'^ + fisx'^)ydx, 
H =\y^ + \-\x^ -tx\ t G (-|,0). 

One can formulate the following conjecture concerning the maximal number of zeroes 
of the generating function I{t) and the corresponding maximal number of limit cycles 
produced by the period annulus (called its cyclicity). 

Conjecture 1. (cf. [H]) The period annulus around the center at the origin in (rO)- 
(r22) has the following cyclicity under small quadratic perturbations: three for cases 
(rl) with a* < a < A, (r3) with | < a < 4, (r4) with A < a <5, (r5) with a = A, (r6) 
with a > 4 and (rlO), and two otherwise. 

We note that a* G (|,3) is determined from a transcendental equation [15j and can 
be calculated numerically, a* = 2.0655... 

Let us mention that for some of the cases Conjecture 1 has already been verified. 
The standard elliptic case (rO) is studied entirely in |T9]. The reversible Hamiltonian 
case (r2) has been investigated in a series of papers, see [1] and the references therein. 
Cases (r5) with 6 = 2 and (rlO) were considered in [H] and [13], respectively. Cases 
(r5) for & 7^ 2, i and (r6) with b G (|, 2) are studied in [2|. Case (rl) with b G (-1, |) 
and b G (|,3) are considered in [22] and [I5], respectively. Below we will concentrate 
our efforts mainly on the cases where Conjecture 1 remains open and will include in 
our lists that follow only these unsolved cases. 

In order to reduce the number of Hamiltonians, we replace 
{x,y) — >• {l/x,y/x^) in cases (r5) and (r6), 
{x,y) (x/(-4t),y/(-4t)3/2) in cases (r7) and (r8), 
(x,y) (2tx, {2tf/^y) in case (r9), 

{x,y,t) {{—6sy^^x,y, (— 48s)"^/^) in cases (rll) and (rl2), 
{x,y) {3tx,9t'^y) in cases (rl3) and (rl4), 
{x,y) — > {1 / {3tx) , y / x^) in cases (rl5) and (rl6), 
{x,y) {~x/3t, —y /3t) in cases (rl7) and (rlS). 
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In this way we obtain a reduced hst of cases to study as follows: 



(rl) 


m 


= Ih=o^ ^ifJ'i + fJ'2X + fj^sx ^)ydx, (6 ^ (-1, 1) U (|,3)) 


(rll) 


m 


= Jfj^Q x{iii + fj,2t^^x^^ + ii^tx^)ydx, {h = ^) 


(rl2) 


I{t) 


= iH=o^ {iJ,i + fi2tx-^ + f^st X )ydx, (0=3) 




TT 

n — 


l.,2 1 3-b 1 6-1™ 1 1-36 ™2 j.™3 
' 2^ 6(6+1) ' 6+1-^ ' 2(6+1)-^ 


(r3) 


m 


= Jh^o x'^^iiii + ii2X^ + ii3X-^)ydx, 


(r4) 


m 


= Jh^o x~'^{Hi + H2X~'^ + iisx^)ydx, 






_ 1 2 1 6+3 1 6-1 2 1 36+1 4 ^3 
- 2^ ^ 24(6+1) ' 4(6+1)-^ ' 8(6+1)-^ 


[r5) 


m 


= Jh^q x{fii + fi2X'^ + fi3x)ydx, {b = |) 



r6) 

r7) 
r8) 



lit 
H - 

I{t 
I{t 



rl3) I{t 

rl4) I{t 

rl5) I{t 

rl6) I{t 



r9) I{t 
rl7) /(t 
rl8) I{t 



= Jh^o{^^i + ^^2X + ^Msx ^)ydx, (6^(i,2)) 

- i„2 I 1-26 2 , 6^™3 I 2-6 4 _ ^ 
" 2^ ^ 2(6+1)-^ ^ 6+1-^ 4(6+1)-^ 

= Jh^o x~^{ij^i + ii2t~^x^ + H3txr ''^)ydx, 
= /^=o ^~^(a*i + l^2tx''^ + iizt~^x^)ydx, 
= JH=oil^^ + A«2^2;~^ + iJizt-'^x^)ydx, 
= /^=o + l^2t~^x'^ + ii^tx~'^)ydx, 

= /^=o + pi'2t~^x'^ + ijLstx~'^)ydx, 

= /^=o ^"^(/^i + l^2tx~^ + iJ,3t~^x'^)ydx, 

l.,2 



JH=oif'i + f'2tx~^ + fj,st~^x^)ydx, 
Ih=o ^^^(a*i + fJ'2t^^x^ + fi3tx^^)ydx, 
Jh=o x^^{fii + ^2tx~^ + fist^^x'^)ydx, 



1^2 -I- i'r3 
2^ -T 



t, te(-io). 



It is useful to know the dimension of the Picard-Fuchs system satisfied by the ba- 
sic integrals Ik{t) — J^^^^ x^ydx involved in the formulas above. Let us rewrite the 
equation H — into the form = T^{x, t). Then for any /c e Z one obtains 



Lit) ^"T^'ydx = x'^ydi^y') = I x'^dy' = -| ^ 



IS{t) 



X 



fe-1 



y^dx. 



Therefore 



X 



k-l 



2k \ 
xV + —V \ ydx 



0. 



Taking V = Ao + A^x + A2x'^ + A^x^ + A4x'^, in terms of the basic integrals Ik{t) the 
last identity is equivalent to 



{2k+12)AJk+3+{2k+9)AsIk+2+{2k+6)A2lk+i+{2k+3)AiIk+2kAoIk-i = 0. (3.3) 



13 



Using (3.3) with A; = 0, 1, 2, . . . and with k — —1, —2, ... we obtain that 

Ik = span(/o, h, h), k> 3; 4 = span(/_i, Iq, Ii, I2), k< -2 {A4 ^ 0). 
Similarly, if A4 = 0, one obtains 

4 = span(/o, 4), k> 2; 4 = span(/_i, 4, 4), k< -2 (A4 = 0, A3 7^ 0), 
and finally, 

4 = span(7o), A; > 1; 4 = span(7_i, Iq), k< -2 (A4 = A3 = 0, A2 7^ 0). 

Here "span" means in general a polynomial -R[t,t~^] span. As a consequence, we can 
formulate a result about the dimension of the Picard-Fuchs system satisfied by the 
basic integrals in (rl), (r3)-(r9) and (rll)-(rl8). 

Proposition 4. In cases (rl), (r3)-(r4) with b = — |, (r6) with b = 2, (r9), (rll)- 
(rl2), (rl7)-(rl8) the Picard-Fuchs system is of dimension 3 while in the remaining 
cases it is of dimension 4- 

To derive the Picard-Fuchs equations, we use (3.3) together with the relations 



We can use (3.4) to further simplify some of the integrals above. Thus we get the 
final fist of reversible cases of genus one yet to study: 




(3.4) 



(rl),(rl2) I{t) = J^^, 
(rll) I{t) = /^^o 





(r3) 
(r4) 




i ( ^ 2_ / 3b+l \ p ('_Q _1 V 

^^1^3(&+1)' 3(6+1) Y 6+3^' 3^' 




te U7rfTT,00 otherwise. 



(r6) 
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(r7),( 


;rl4) 


I{t) 


= Ih 


(r8) 




I{t) 


= Ih 


(rl3) 




I{t) 


r 

= Jh 






T(f) 


- f 

~ Jh 




'r1 fi) 


T(f) 


- f 
~ Jh 
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2y 


(r9) 




m 


= Ih 


(rl7) 




i{t) 


= Ih 


(rl8) 




i{t) 


= Ih 






H = 





j_ 

12- 



2'" '3 ' " ^ V 6' 

We point out that the above formulas are related in each case to the period annulus 
around the center at (1, 0) obtained for t e {tc, ts) where tc is the level corresponding 
to the center and ts is the level of the contour at which the annulus terminates. 

4. The generating function in the Lotka-Volterra case 

As in the previous section, we will obtain the complete list of generating functions 
I{t) for the generalized Lotka-Volterra systems (2.1) in the cases when all their orbits 
are elliptic curves. Apart of the reversible case, the equivalence classes for the Lotka- 
Volterra case listed in Propositions 1 and 2 are a result of affine transformations. 
Therefore we can choose by one representative from each such a class and then use 
the generating function corresponding to it. We will always choose the cases whose 
parameters satisfy conditions (I) and (II) above. Take real coordinates z = x + iy 
and consider a small quadratic perturbation as in (3.1). In the reversible case (2.4), 
the first integral and integrating factor are respectively [1] 



A(A-1)V,, A + 2-, , „ 




i/(.Y.,) = .r|^^-^^^A-- — ) I. M = .Y-' (4.1) 

where 

h -\- I 1 4r 

A = r^, 6^0,±l,-, X = l- 



6-1' ' ' '3' A-1 
The generating function for small quadratic perturbations reads [H], Theorem 3: 

I{t)= I x^'^[^iy + H2yx~^ + ^J^^{x - l)y'^]dx. (4.2) 

JH=t 

In the Hamiltonian triangle case (rlvl), the integral I{t) takes another form and this 
case has already been studied in p^. For this reason (rlvl) is not included in our 
list below. By (4.1), (4.2), (2.6) and the formulas of first integrals of genus one 
(rlv2)-(rlv6), we get the list of explicit expressions of I{t) for all other reversible 
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Lotka-Volterra cases as follows: 



(rlv2) 


m 


lH=t 




l)y ^)dx, 




H = 




t6(-|,0) 




(rlv3) 


lit) 




[fiiy + fi2yx ^ + fi3{x-l)y ^)dx, 




H = 


x'ib' - 


-ix-2f), te{-l,0) 




(rlv4) 


m 


~ lH=t ^ 


~^ifiiy + fi2yx~^ + fisix - 


l)y ^)dx, 




H = 




-ix-l)'), tG(-i,0) 




(rlv5) 


lit) 


= lH=t ^ 


"^(/iiy + fl2yX~^ + fl3ix - 


l)y^^)dx 




H = 


— - / 1 2 
X 2(ir 


+ + tG(f,oo) 




(rlv6) 


m 


= lH=t ^ 


"^(/iiy + H2yx~^ + fisix - 


l)y~^)dx 




H = 


— i / 1 2 

X 2[-y^ 


+ (x + 3)2), t G (16, oo). 





Clearly the elliptic curves in cases (rlv5) and (rlv6) are obtained by introducing a new 
variable x — >■ so that they are given by the level sets of H{x'^, y). Let us also recall 
that the functions I{t) in (rlv2)-(rlv6) are the coefficients at in the expansion of 
the displacement map obtained for the special perturbations in (3.1) which keep the 
coefficient at e zero. 

What concerns the small quadratic perturbations of the generic (nonreversible) Lotka- 
Volterra system, in (3.1) we have M = x^~^y^^^ (modules needed outside the first 
quadrant but we will omit them) where H is determined by (2.2). In the generic case 
I{t) is given by the complete elliptic integral [14j, Theorem 2 

I{t)= I I x^-^y'''\ni + fi2X~^ + fi3y~^)dxdy. (4.3) 

J Jlnt5{t) 

Using (4.3) and the list of first integrals of genus one (lvl)-(lv5), we easily obtain the 
concrete form of I{t) for all the CcLSGS cLS follows: 



(Ivl) 


I{t) 


= Ih- 




ly 


l(/il +/i2X ^ + fi^y ^)dx, 




H = 


2 1 

X3y3\ 


(1- 


- X ■ 


-y), tG (0,432-1/3), 


(lv2) 


m 


= Ih- 


-t^ 


7 

ey 


"^(/ii + /i2X~i + H3y~^)dx. 




H = 


1 

X ey 


-( 




x-y), te (4321/6,00), 


(lv3) 


m 


= Ih- 




7 

ey 


"^(/ii + fi2X~^ + fi3y~^)dx. 




H = 


x'ly 


-H 




x-y), te (432^/^00), 


(lv4) 


m 


= Ih- 


-t^ 


-iy 


"^(/ii + /i2X-i + H3y~^)dx. 




H = 


1 

X -iy 


-H 




x-y), te (2v/2,oo), 


(lv5) 


lit) 


= Ih- 


■-t ^ 


4 

:iy 


"^(/ii + /i2X~i + fi3y~^)dx. 




H = 


1 

X 3y 


-H 




x-y), te (432i/^cx)). 
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As in Section 3, we formulate a general conjecture about the cyclicity of the period 
annulus in the Lotka-Volterra systems having all their orbits elliptic or conic curves. 

Conjecture 2. The cyclicity under small quadratic perturbations of the period an- 
nulus surrounding the center at the origin in the Lotka-Volterra system (2.1) is as 
follows: three in case (rlvl) {the Hamiltonian triangle) and two in all other cases 
(rlv2)-(rlv6), (lvl)-(lv5). 

Except for the isochronous center Si and the Hamiltonian triangle, it is very likely 
that Conjecture 2 is still open. We recall that quadratic perturbations of the general 
Lotka-Volterra system have been considered in [23]. However, in the recent book [1], 
page 379, V.I. Arnold declared that the problem with the Lotka-Volterra system still 
remains open. 

Remark about the codimension four case. In this remark we discuss in brief 
the generic (non-reversible) codimension 4 center. In complex coordinate z = x + iy, 
the related system becomes 

z = -iz + 4:Z^ + 2\z\'^ + az^ , a G C \ M, \a\ = 2. 

In {x,y) = (X^,y) coordinates, where 

4 

Y = (2 + b)x + cy, X = 1 + 8x H -Y^, a = b + ic, 

V J 2 + 5 

the system has a first integral of the form [Hj (below the bars are omitted) 

It is seen that the level sets of this first integral are (generically) elliptic curves. 
Therefore the generating function I{t) whose zeroes correspond to limit cycles in 
the perturbed system are given by the following complete elliptic integral (cf. [H], 
Theorem 2 (iii)) 

H^) = I I x~^[fii + fi2y + fisV^ + fJ'ii'^^y'^ - x^)]dxdy, k = 

J JH{x,y)<t Z + 

The conjecture that the integral J(t) has at most three zeroes in the interval (tc, tg) = 
(— i2(2-b) ' ^) corresponding to the period annulus around the center at (xc, yc) = (1, 0) 
is still open. 

5. Zeros of Abelian integrals for some of the reversible cases 

In this section we study the generating functions I(t) related to the reversible quadratic 
systems (rl8) and (rll). These systems contain no parameter, have a unique period 
annulus around the center at (1, 0) and the Picard-Fuchs system satisfied by the com- 
ponents of J{t) = I'{t) is of dimension three and two (respectively). Our main result 
is the following. 
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Theorem 3. The exact upper bound of the number of the limit cycles produced by the 
period annulus under quadratic perturbations of the reversible system (rl8) or (rll) 
is two. 

To prove Theorem 3 we study first one-parameter analytic quadratic perturbations 
of (rlS) and (rll). According to the formulas derived in Section 3, the number of 
the limit cycles of such perturbations is bounded by the number of the zeros of the 
following generating functions 



I{t)^ j^^^{lii + li2X-' + ii3x)ydx, te(^-^,Q^, Ht^^y^-^x^ + ^x^-t, (5.1) 
for system (rl8) and 

I{t) = J^^^x-\iii+fi2X-^+fi3x)ydx, t e (^-^'O) ' = ^y^+^-^x-tx\ (5.2) 

for system (rll). In the above integrals S{t) denotes the unique oval of the real 
algebraic curve {{x, y) e : Ht(x, y) — 0} for a given t e (— |, 0) and jii are arbitrary 
real constants. The ovals 5{t) in (5.1) form a period annulus which is bounded by the 
homoclinic loop {Hq = 0} going through the saddle at the origin, while in the case 
(5.2) they form a period annulus bounded by the parabola {Hq = 0}. In both cases, 
the ovals S{t) exist for t e (— 1,0), I{t) is analytic in a neighborhood oi t — —1/6, 
and /(— 1/6) = 0. Consider the derivative J{t) — I'{t), see formulae (5.12), (5.13) 
bellow. The key ingredient of the proof of Theorem 3 is the following result. 

Theorem 4. The three-dimensional vector space of Abelian integrals J{t) = 

t e [— |, 0), defined by (5.1) or (5.2) is Chebyshev. This means that each integral J{t) 

has at most two zeros {counted with multiplicity) in the interval [— |, 0). 

Before proving Theorem 4 we need some preparation. 

5.1. Monodromy of the level curves. Let us denote 

Vt = {{x,y) e C^ Hi{x,y) = 0}, Vt ^ {[x : y : z] e CP^ Ht (-, ^) = 0}. 

If t ^ -|, 0, then both Vt and are smooth curves and Vt is a compact Riemann 
surface of genus one which is also the compactification of the affine elliptic curve Vt- 
We have = U oo where oo = [0 : 1 : 0] and therefore 

rank Hi{Vt, Z) = 2, rank Hii^t, Z) = 2. 

For (5.2), all compact complex level curves Ft have two complex conjugate common 
points, namely: P+ = {Q,i^J^) and P~ = (0, — iy^) in affine coordinates, as well 

as one common point at infinity. For (5.1) there is one common point, the point 
at infinity. For both cases, when blowing-up these points on the complex projective 
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plane P^, one obtains a compact smooth rational surface S and an analytic map 
5* — >■ P where the projection vr is induced by one of the rational maps 




(5.3) 



2^^ +3 Ix 



(5.4) 



In both cases 5* is an elliptic surface in the sense of Kodaira [16] with three singular 
fibers 7r~^(0), 7r~^(— 1/6), 7r~^(oo). In particular, this allows one to compute the global 
monodromy of the related homology bundle as described below. 

5.1.1. The monodromy of the fibration (5.3). Denote by 0^ the two points 
on with coordinates {x = 0,y = ±\/2i) and let Tt = Tt\ {0^,0"}. We shall 
determine the monodromy of the homology bundle associated to (5.3) with fibers 
Hi{Tt, Z) = Z^. Let S{t), 7(t) be a continuous family of cycles, vanishing at t = — 1/6 
and t = respectively, and let a(t) be a cycle represented by a small loop around 0"*" 
on Ff. We need to describe the monodromy of these cycles on the plane C\{0, —1/6}. 

The precise definition of the families S{t), jit) is as follows. For t G (—1/6, 0) the 
polynomial — — t has three real roots xi(t) < < X2(t) < xs(t). Let / be a 
simple loop on C \ {xi{t) , 0, X2(t) , x^it)} which makes one turn about xi(t), 0, X2(t) 
in a positive direction and does not contain in its interior x^lt). 

We define •yit) to be the cycle represented by / on F^ (that is to say by one of the 
two pre-images of / under the projection (x, y) ^— x). The cycle ■y(t) is defined up to a 
sign. The cycle ±5(t) is defined in a similar way, but the simple loop / is supposed to 
make one turn about X2{t),X3{t) in a positive direction, and does not contain Xi{t), 
in its interior. 

We note that xi(0) = X2(0) = 0, ^3(0) > and hence •yit) is a cycle vanishing at 
t = 0. Similarly, 6{t) is a vanishing cycle at t = —1/6, as xi(— 1/6) < 0, X2(— 1/6) = 



be the monodromy representation related to the elliptic fibration associated to (5.3). 
The image 



is a group generated by 1^,11,11^ = II o /g, the monodromy operators corresponding 
respectively to the oriented loops /o, ^1,^00 = h° ^o- Here /q is a simple loop which 
makes one turn about in a positive direction and does not contain —1/6 in its 
interior. Similarly li is a simple loop which makes one turn about —1/6 in a positive 
direction and does not contain in its interior. 

Lemma 1. In the global basis {a{t),6{t),'~f{t)) {with appropriate orientation of the 



X3(-l/6) >0. 



Let 



n : 7ri(Cr \ A, to) ^ Aut {Hi{Tt, Z)), A = {-|, 0, 00}, to e (-|, 0) 



n(7ri(Cpi\A,to)) 
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cycles), the monodromy operators associated to the fibration defined by (5.3) are: 



1*0 -Mo 



/-I 1 o\ 
1 



V 



-1 1 



/* : Ml 



1 
1 1 
1 



r : Moo = MiMo 



/-I 1 l\ 
1 
-1 1 



Proof. The identity 



2m — = — 



dx 1 

.(t) xy 

shows that lQ{a{t)) = —a{t). The remaining claims follow from the usual Picard- 
Lefschetz formula (it is enough to describe the monodromy of the roots Xj(t)) and 
therefore the details are omitted. □ 

5.1.2. The monodromy of the fibration (5.4). Here we determine the mon- 
odromy representation 

n : 7ri(CP^ \ A, to) ^ Aut {Hi{Tt, Z)), A = {-|, 0, oo}, to e (-|, 0) 

associated to the elliptic fibration defined by (5.4). The automorphisms IqJU^^ 
associated to the oriented loops lo,h, loo are defined as in the preceding subsection. 

For (5.4) the global monodromy does not follow from the Picard-Lefschetz formula. 
The local monodromy (around the singular fibers), however, depends only on the 
topological type of these fibers and is computed in [16j. The topological type of the 
fibers 7r~^(0), 7r~^(— 1/6), 7r~^(cx)) on its turn is computed in [S] and it is respectively 
(III), (Ii), (IV*) (using the Kodaira notations, see e.g. [171 Table 6.1]). We conclude 
that, up to a conjugation, the monodromy operators are given by (e.g. [TOl Table 1]) 



/* • M 



II : Ml 



1*0 -Mo 




By abuse of notation, we denote by 6(t) G Hi(Tt,Z), t G (—1/6,0), the continuous 
family of cycles, vanishing when t — — | along a path connecting t to — | in (— |, 0). 
Let 7(t) = lQ{5{t)). The continuous families of cycles 5(t),7(t) are defined in C \ 
{-i,0}. According to [IH], (5(t),7(it)) is a basis of H^(Tt,li). 



20 



Lemma 2. In the global basis {5{t),'y{t)) {with appropriate orientations of the cycles), 
the monodromy operators associated to the fibration defined by (5.4) are: 



Proof. For the global sections 5,7 of the homology bundle we have l^S = 7, llS = S 
and hence 

/ * \ f 1 * 

Mo = , Ml = 

\^ 1 * y \ * 

The relations Tr{Mi) = 2, Tr(Mo) = 0, det{Mo) = 1, Tr(MiMo) = Tr{M^ 
imply the result. □ 



5.2. Wronskians. Let 

x''dx 
ujk = , k E Z 

y 

be polynomial one-forms on C^. They induce meromorphic differential one-forms on 
the compact Riemann surface which are denoted in the same way. 
Let us introduce the Wronskian 



Js{t) Jj{t) J5{t) J-y{t) 

where, by abuse of notation, S{t),''y{t) are continuous families of closed loops, with 
the properties : 

• they intersect transversally at a single point 

• Fj \ {(5(t),7(t)} is homeomorphic to a rectangle whose opposite sides are iden- 
tified to 6(t) and ■y(t) 

• The poles of Uk are contained in the interior of this rectangle 

As LJo is holomorphic in both cases, the above Wronskian can be computed by making 
use of the reciprocity law for abelian diff'erentials of the second and third kind ^ p. 
647] as shown next. Below, we denote by C a certain nonzero constant. In the case 
(5.2) the Wronskian is a rational function and 

Ws,^{u;k, ujo) = = <^ 0, k = 2, (5.5) 



c 
t ' 


k = 


1, 


0, 


k = 


2, 


c 

i2' 


k = 


3. 
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In a similar way, in the case (5.1) we have 



= a,„ + j^^ = ^ (5.6) 

and 

Ws,,{uj^,uj,) = ^^^^^ = C ^ 0. (5.7) 

5.3. Asymptotics of the Abelian integrals. Here we study, for suitable /c, the 
asymptotical behavior of Jkif) = /^^j) (-^fe, where S{t) is the continuous family of cycles 
vanishing at the center (associated to t — —1/6), near t — oo and t — 0. Below, 
we shall write J{t) {t — toY log(i — ^o)^, r, s e M provided that for every sector 5" 
centered at the critical value & A there exists a nonzero constant Cs such that 
\J{t)\<Cs\t-to\'\logit-toY\. 

5.3.1. The case (5.1). Here we have 

^ , , /' x'^dx f x^dx , „ , 

Jk{t)- / = / ; =, A; = 0,±1. 

-Im y -Im J-lx^ + a;2 + 2t 

(a) Near ^==00, the change of variables x — t^u leads to 

^ , , k 1 f u'^du k 1 ~ , , 

When \t\ 00 (with bounded argument) the integral Jk{t) tends to a finite 
constant. Consequently 

Jk{t)<t~'~K keZ. (5.8) 

(b) Near t — 0, the Abelian integral Jk{t) can be expanded as follows 

Jk{t) = / ^k-l f + Q{t) 

^TTZ Jy(^t) ^ Ja{t) 

where Q is a meromorphic function in a neighborhood of t = (this follows 
from Lemma 1). Therefore 

J,{t)^^^+P{t)log{t) + Q{t) 

where P, Q are meromorphic in a neighborhood of t = 0. We claim that 

log(i) if A; = 0, 
Jk{t) <{ 1 if A; = 1, (5.9) 
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Indeed, for A; > 1 

lim Jk(t) = 2 / 
^0- Jo 



where ^3(0) = 3/2. For A; = the integral Jo{t) is of the first kind and its 
behavior is well known. For k — —1 the leading term of the expansion of J-i{t) 
is given by 

dx 1 dx 



ly^ x\/ x'^ — t \/i Ji x\/x^ — 1 
(c) Near t — —1/6 the functions Jk{t) are holomorphic and Jk{t) ^ 1. 

5.3.2. The case (5.2). We have 

, x'^dx f x^dx 
Jk{t)= =/ , A; = l,2,3. 



r x'^dx f 

Jm y Js{t) 



+ X-1 



(a) Near t — 00, the change of variables x — t 3U leads to 

Jk{t) = t-— / , = = t-—Jk{t). 

Js{t) ^2u^ + t--3u-l 

When |t| ^ 00 (with a bounded argument) the integral Jfc(t) tends to a finite 
constant (possibly zero) and hence 

4(^)<r^, kez. (5.10) 

(b) Near t = any Abelian integral of the first or second kind is a power series in 

t^^^ (because the eigenvalues of 1^ are fourth roots of the unity). As t tends to 
two roots of te^ +a:/2 — 1/3 tend to 00 which suggests us to make the change 
of variables x — f^u. Then we have 

Jk{t) = r(2+4) / = ^(2+4) j,(t) 

Js(t) ^2u^ + M - |t5 

and hence 

Jk{t)<t-^-^^'-\ kez. (5.11) 

(c) Near t — —1/6 the integrals Jk{t) are holomorphic and Jk{t) < 1. 



5.4. Proof of Theorem 4. We have to prove that the derivative 

c'^dx f xdx f x^c 

h/i2 / h/i3 / 

's(t) y Jsit) y Js{t) y 
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J{t)^I'{t)^^iJ ^ + ^iJ ^+nJ ^ (5.12) 

J5(t) y Js(t) y Js(t) y 



(in the case (rll)), or 



J{t) = I'{t)=^^J - + - + /^3/ ^ (5.13) 

J Sit) y Jsit) xy J Sit) y 

(in the case (rl8)) has at most two zeros. For this, we use the method of Petrov, see 
e.g. [20], based on the argument principle. 
Introduce the function 

Jow Js{t) y 

It is real analytic on (—1/6, 0) and has an analytic continuation in the complex domain 
V = C\[0, oo), because the period Jo(t) of the elliptic curve ft does not vanish there, 
including at the point t = — 1/6. To bound the number of the zeros of F in D it is 
enough to find this number in the smaller domain = T) Cl {t : r < \t\ < R} for 
r sufficiently small and R sufficiently big. We are going to evaluate the increment of 
the argument of F along the boundary of "P/jr, oriented in a positive direction. 

5.4.1. Zeros of F in the case (rl8) and (5.1). Along the boundary of the 
small circle {|t| = r}, according to (5.9), the increase of the argument of F is at 
worst close to vr, and along the boundary of the big circle {\t\ = R} this increase 
is at worst close to 27r/3, see (5.8). Denote by F"^ the restriction on (0, cx)) of the 
analytic function obtained as an analytic continuation of F along an arc contained in 
the upper (lower) complex half-plane Im ± t > respectively. The function F is real 
analytic on (— oo, 0) which implies that along the interval (0, oo) the imaginary part 
Im(F+(t)) of F+ satisfies 



2i Im (F+ (t) ) = F+ (t) - F+ (t) = F+ (t) - F~ (t) . 

Assume for a moment that F~^{t) —F~{t) has at most one simple zero on the interval 
(0, oo). Then summing up the above information we conclude that F{t) has at most 
two zeros in the complex domain T>jir, and hence in T>. This result is obviously exact. 

It remains to prove that F~^{t) — F~{t) has at most one zero on (0, oo). Clearly 
F~ is an analytic continuation of F^ along an arc making one turn about t = 
in a positive direction. This shows that F^{t) — F~{t) is obtained as an analytic 
continuation of the function 

F{t)-ll{F{t)) = ^- = ^mmf^^^o) ^ tG (-1/6,0) 

h{t) ^0 Jl*5it) ^0 Jsit) ^0 Ji*s(t) ^0 

along an arc contained in the upper complex half-plane. Here u = fj,iU!Q+^2^^-i+fJ'3^i, 
/o5(t), 6{t) are cycles on the elliptic curve defined by Ht = 0, with two removed points 
0^, and l*Q5{t) = 5{t) - -f{t) + a{t). Let 

iW{t), tG(0,oo) 
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be the real analytic function obtained as an analytic continuation of 



iWs(t),i*s{t){^^,^^o), te (-1/6,0) 
along an arc contained in the upper complex half-plane. As 

Wit) 



F^t)-F-{t) 



Is{t) ^ol 



then we shall show that iW{t) has at most one zero on (0, oo). We use once again 
the argument principle. We shall show that the analytic continuation of W{t) in 
V — C \ (—00, 0) has at most one zero counted with multiplicity. For this purpose we 
consider the complex domain 

VRr = Vf\{t: \t\ > r, \t\ < R}. 

As before, denote by W^{t) the restriction on (— cxd,0) of the analytic function ob- 
tained as an analytic continuation of W{t),t G (0, 00) along an arc contained in the 
upper (lower) complex half-plane Im ±t > respectively. Along the interval (—1/6, 0) 
we have 

W+{t) - W-{t) = Ws(^t),^mi^^^o) - {l*o)-'Wsit),roSit)i^,^o) 

= W5(t),l*Sit) (i^, t^o) - W(l*)-is(t),Sit) (<^, <^o) 

= W^<5(i),5(t)+a{t)-7(t) 1^0} - Ws{t)+a{t)+'r{t),d{t) {^^, ^o) 

= 2iyj(t),„(t)(a;,a;o) 



UJo- 



2t Js{t) 

Note that, according to section 5.1.1., the functions W^{t) are single- valued in a 
neighborhood oi t = —1/6. Therefore 



W+{t)-W-{t) = -^ [ uo 
V 2r J5{t) 



2t J5{t) 

along the interval (—00, 0). We conclude that the imaginary part i{W^ {t)—W~ {t)),t e 
(—00, 0) of the analytic continuation of iW{t),t e (0, 00) does not vanish. 

It follows from the asymptotic expansion of near t — Q obtained in 5.3.1.(b) 
that 

Therefore along the small circle {\t\ — r}, the increase of the argument of W{t) is 
at worst close to vr. Along the border of the big circle the decrease of the argument 
is close to 27r/6, sec (5.8). Summing up the above information we conclude that the 
total increase of the argument of W{t) along the border of X^^^ is close to (or smaller 
than) 47r — 27r/6 and hence W{t) has at most one zero. Therefore the maximal number 
of the zeros of F{t) in T> is two [and this result is exact). 
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5.4.2. Zeros of F in the case (rll) and (5.2). In the same way as above we may 
study F(t) in the case (rll). This leads, however, to a bound of the number of the 
zeros in [—1/6, 0) equal to three. To improve the bound we consider first the function 
F{t) = F(t) + /i4, where 114 is a real constant. We shall show that F{t) has at most 
three zeros in the complex domain T> = <C\[0, 00). For this we consider once again 
the domain T>Fir = T> (1 {t : r < \t\ < R} and evaluate the increase of the argument 
of F along its border. Along the boundary of the small circle {\t\ = r}, according to 
(5.11), the increase of the argument of F{t) is close to Sir, and along the boundary 
of the big circle {\t\ = R} this increase is close to 0, see (5.10). Along the interval 
(0, 00) the imaginary part of F is equal to 

Wit) 



where, according to (5.5), for suitable constants ci,C2, 

Cit + C2 



W{t) 



It follows that the imaginary part of F{t) has at most one zero along (0, 00). The 
argument principle implies that F{t) has at most three zeros in the complex domain 
V. Suppose now that for some /ii,/i2,At3 the function F(t) has exactly three zeros in 
V, and hence in the domain V^ir for r sufficiently small and R big enough. According 
to (5.10) for |t| sufficiently big we have 

F{t) = cr''/^ + o{\t-''/^\). 

It follows that for sufficiently small /i4 a real zero of F{t) bifurcates from 00 on the 
interval (— oo,0) and this zero is not contained in the domain Therefore F{t) 

will have at least four zeros in the domain V, in contradiction with the result proved 
above. Therefore the maximal number of the zeros of F{t) in T> is two {and this result 
is exact). 

Theorem 4 is proved. □ 

5.5. Proof of Theorem 3. Denote the open period annulus of the fixed reversible 
system (rl8) or (rll) by 11. Let X\, A G A be the set of all quadratic plane vector 
fields, analytic with respect to A and such that Xq coincides with (rl8) (respectively, 
with (rll)). Theorem 3 can be reformulated as follows: 

The cyclicity Cycl{Il, Xx) of the open period annulus U is equal to two 



which means that Xx has at most two limit cycles which tend to 11 as A tends to zero. 
It is shown in [91 Theorem 1] that if the cyclicity Cycl{Il, Xx) is finite, then there 
exists a germ of an analytic curve A(e), such that 

Cycl{U,Xx) = Cycl{U,Xx{6)). 
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In other words, it is enough to study one-parameter deformations. Let I(t) be the 
first non-zero Poincare-Pontryagin (or generating) function associated to such a per- 
turbation. The derivative J(t) = /'(t) is of the form (5.12) or (5.13) and we proved 
that this function has at most two zeros (Theorem 4), and hence I{t) has at most 
three zeros on [—1/6,0), one of them being t = —1/6 . Therefore, by a standard 
argument, the cychcity of the open period annulus of the perturbed one-parameter 
quadratic system is at most two. 

It remains to show that the cyclicity Cycl{Il, Xx) is finite. It is shown in [9j that if 
Cycl(Il, Xx) = oo, then there exists a Poincare-Pontryagin function I{t) (associated 
to some one-parameter deformation) with infinite number of zeros. As I{t) is an 
Abelian integral, then this is clearly impossible. This completes the proof of Theorem 
3. □ 
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